Can stochastic quantization evade the sign problem? 
— the relativistic Bose gas at finite chemical potential 
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A nonperturbative study of field theories with a complex action, such as QCD at finite baryon 
density, is difficult due to the sign problem. We show that the relativistic Bose gas at finite chemical 
potential has a sign and 'Silver Blaze' problem, similar to QCD. We then apply stochastic quantiza- 
tion and complex Langevin dynamics to study this theory with nonperturbative lattice simulations. 
Independence of chemical potential at small and a transition to a condensed phase at large chemical 
potential are found. Lattices of size iV 4 , with N = 4,6,8,10, are used. We show that the sign 
problem is severe, however, we find that it has no negative effect using this approach. This improves 
the prospects of applying stochastic quantization to QCD at nonzero density. 
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Introduction - Field theories with a complex action are 
difficult to treat nonperturbatively. Because the weight 
■ in the path integral e~ s = \e~ s \e l<p is not real, standard 
7— I \ numerical approaches based on a probability interpreta- 
tion and importance sampling cannot be applied. This 
has hindered progress in condensed matter and many- 
body theories, such as frustrated quantum spin systems 
and strongly correlated electronic systems away from 
half filling. An outstanding example is QCD at nonzero 
baryon density. Several methods have been devised to 
circumvent the sign problem in QCD: reweighting 0,0], 
Taylor series expansion [^, 0| , imaginary chemical poten- 
tial and analytic continuation Q , and the use of the 
canonical ensemble 0, H| and the density of states [§] . In 
general these methods can only be applied in a limited re- 
gion of the phase diagram (such as high temperature and 
small chemical potential), suffer from overlap problems, 
and/or are restricted to small volumes. 

At vanishing temperature, so far all methods have 
failed. It is well-known what is supposed to happen: if 
\x is the quark chemical potential (so that the fermion 
determinant satisfies detM(/i) = [det M{— //)]*), one ex- 
pects a transition to a condensed phase (nuclear matter) 
at a critical chemical potential fi c ~ mjv/3, where win is 
the nucleon mass. Below fi c , physical observables should 
be completely independent of \x (at strictly zero temper- 
ature) , even though microscopically /i explicitly enters in 
the Boltzmann weight. This exact cancelation is known 
as the 'Silver Blaze' problem [lOj and is intimately tied to 
the sign problem, as has been demonstrated in Random 
Matrix Theory [n], EH ■ It is an outstanding challenge 
to probe the QCD phase diagram in the low-temperature 
region at nonzero chemical potential. 

The sign problem in QCD arises not because of the 
anticommuting nature of fcrmions. Instead, it is due to 
the behaviour of the action under complex conjugation 
when the chemical potential is nonzero. Therefore, this 
sign problem can also be studied in bosonic field theories 
with a chemical potential coupled to a conserved charge 



and an action satisfying S(fi) = [S(—fJ>)]*. These theories 
suffer from exactly the same sign and Silver Blaze prob- 
lems as QCD. In this Letter we consider the relativistic 
Bose gas at finite chemical potential and study its Sil- 
ver Blaze problem with nonperturbative lattice simula- 
tions employing stochastic quantization 
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Stochastic 

quantization does not rely on the interpretation of e~ s 
as a probability weight. Instead the proper distribution 
is obtained as the equilibrium distribution of a stochas- 
tic process, described by a Langevin equation. In the 
case that the action is complex, the Langevin dynamics 
is complexified llotllol . See Ref. 17 1 for a comprehensive 
review and Ref. [l8| for an early application to the SU(3) 
spin model at finite chemical potential. In the past com- 
plex Langevin dynamics has been hindered by numeri- 
cal instabilities and uncertainty about convergence, see 
e.g. Ref. [13]. Recently some of these problems were al- 
leviated by the use of more refined L ang evin algorithms 

23, EL 



22 



231 ] . In particular, in Ref. [23j the method has 
been applied to several models at finite chemical poten- 
tial. In U(l) and SU(3) one link models, where the sign 
problem is mild, excellent results have been obtained. 
For QCD with static quarks, first results on a 4 4 lattice 
indicate a transition from a low-density confining to a 
high-density deconfining phase. Even though the sign 
problem appears to be severe, observables are under con- 
trol. In this Letter we provide considerable evidence that 
stochastic quantization evades the sign problem and is ca- 
pable of handling the Silver Blaze problem, in the context 
of the relativistic Bose gas at finite chemical potential. 

Bose gas at finite chemical potential - We consider a 
self-interacting complex scalar field in d = 4 euclidean 
dimensions, with the lattice action 



(2d + m 2 ) <f>*J x + \{(t>l(t> x f 
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The lattice spacing a = 1 and we take m 2 > 0. The 
lattice four- volume is O = N^N T , where N s (N T ) are the 
number of sites in a spatial (temporal) direction. As 
usual, chemical potential /i is introduced as an imag- 
inary constant vector potential in the temporal direc- 
tion. Note that S(fi) = [S(—fi)]*. At zero tempera- 
ture and in the thermodynamic limit, bulk physical ob- 
servables are strictly independent of the chemical po- 
tential as long as /Li < /x c , with \x c the critical chemi- 
cal potential. At \i = /i c , one expects a second order 
phase transition to the Bose condensed phase, where 
the density (n) = (1/0)9 In Z/d/j, is nonzero. Ignoring 
interactions, the critical chemical potential is given by 
= 2arcsinh(m/2) (corresponding to = m in 
the formal continuum limit). Interactions are expected 
to increase this value [29(. The exact /i- independence 
when fj, < He, even though microscopically fi is mani- 
festly present, is the Silver Blaze problem. 

Complex Langevin dynamics - To apply stochastic 
quantization, we start with the Langevin equation 
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5S\, 



(2) 



where 9 is the Langevin time and r\ is Gaussian random 
noise. Since the action is complex, stochastic quantiza- 
tion relies on a complexification of the fields. The original 
field is first written in terms of two real fields <p a (a = 1, 2) 
(0i + z0 2 ) • In terms of these the action reads 
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Here the antisymmetric tensor e a b, with £12 = —£21 = 
1, £n = £22 = 0, is introduced and summation over 
repeated indices is implied. The term proportional to 
sinh fj, causes the action to be complex. 
These real fields are now complexified as 

0„-^ + i# (a =1,2), (4) 

and the complex Langevin equations read 

O + l) = + (5) 

< x (n+l) = <» + eO). (6) 

Here Langevin time is discretized as 9 = rie, with e 
the time step, the noise rj is real and Gaussian, with 
(r]a,x(n)) = 0, (ria,x(n)i] b ^(n')) = 25 nn >8 ab 8 xx > , and the 
drift terms are determined by 
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FIG. 1: Real part of the density, Re (n), as a function of chem- 
ical potential for lattices of size N 4 , with N = 4, 6, 8, 10. The 
parameters are m = A = 1 and stepsize e = 5 x 10~ 5 . The 
inset shows a blowup around the transition. In the thermody- 
namic limit, the density vanishes below the critical chemical 
potential. 



Explicitly, they read 
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Observables are written in terms of the complexified 
fields using the replacement (|1]). For instance, the square 
of the field modulus is given by 



~ 2 K 



(11) 



and the density is given by n — ^ Ylx Ux > w i^h 
n x = (S a b sinh \i - ie ab cosh /1) 4> a ,x<f> btX+ 4 (12) 
(S ab sinh /i - ie ab cosh /i) (^x^ x+ 4 " ^x&l^+i 

(13) 



' 'a, x4 'b,x+4 



All observables now have a real and imaginary part. 

Results - We have solved Eqs. ([5j numerically, us- 
ing a Langevin stepsize e = 5 x 10 -5 , a total number of 
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FIG. 2: Real part of the square of the field modulus, Re (|0| 2 ), 
as a function of /i. The inset shows a blowup at smaller fi. 

5 x 10 6 Langevin timesteps, and m = A = 1. Lattices 
of size N 4 , with N = 4,6,8,10 were used. No instabil- 
ities or runaway solutions were encountered. The real 
part of the density is shown in Fig. [T] for chemical po- 
tentials < ji < 1.7. Imaginary parts of observables 
were found to be consistent with zero, as they should. A 
transition between a zero-density phase and a condensed 
phase with nonzero density is clearly visible. The inset 
shows a blowup of the transition region. Nonanalyticity 
associated with a phase transition can only occur in the 
thermodynamic limit. We observe that with increasing 
four-volume the transition becomes sharper and the den- 
sity goes to zero below [i c 1.15, as is expected in a 
second order phase transition [3(J. Finite size effects are 
comparable to what is found analytically in the noninter- 
acting system, when < \yP c \. In Fig.[2]we show the real 
part of the square of the field modulus (fTT]l . Again we 
observe /i-independence below \i c and a sharp rise above. 
We note that the value at = is obtained using real 
Langevin dynamics and is therefore theoretically well es- 
tablished. We conclude that the Silver Blaze problem 
poses no obstacle for stochastic quantization and that 
there is no problem in taking the thermodynamic limit. 

Some insight into why this method works can be ob- 
tained by ignoring the interactions (A = 0) [25j| . In this 
case the complex Langevin equations can be solved an- 
alytically and convergence to the exact results can be 
proven, provided that \fi\ < (= 0.9624 for the pa- 
rameters used here). For larger /j,, the free theory is un- 
stable. As we have shown here, interactions shift the 
critical chemical potential and remove the instability. 

Sign problem and phase quenching - To quantify the 
sign problem, we write e~ s = e~ SR ~ lSl = \e~ s \e ltp . In 
reweighting, the phase factor e lv is combined with the 
observable while simulations arc performed in the phase- 
quenched theory, obtained by ignoring the imaginary 
part of the action. Expectation values in the full the- 




FIG. 3: Density in the phase-quenched theory, (n) pq , as a 
function of fi. The inset shows a blowup at smaller \i. 
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FIG. 4: Average phase factor in the phase-quenched theory, 
Re (e^) pq , as a function of u, indicating the severeness of the 
sign problem in the thermodynamic limit. 

ory are then reconstructed as (0)f u ii = (Oe lv ) pq / (e iip ) pq , 
where the subscript pq indicates phase quenched. The 
average phase factor (e lv ) pq is the ratio of two partition 
functions, ( elip ) pq = Zf u u/Z pq = e~ nA ? , and vanishes 
exponentially in the thermodynamic limit CI — > oo, when 
0. Here A/ is the difference between the free energy 
densities in the full and phase quenched theories. 

We have studied the phase-quenched theory numeri- 
cally. In this case, real Langevin dynamics can be em- 
ployed, using only Eq. ((SJ) with $ a = 0. In the expression 
for the density (|12p , only the term proportional to sinh fi 
is preserved. Simulations were carried out using the same 
parameters as above. No difference in convergence prop- 
erties were found. The density in the phase-quenched 
theory is shown in Fig.[3J We observe that the density in- 
creases linearly as soon as /x ^ 0. This is not unexpected, 
since the density is proportional to sinh/i. However, it 
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makes it even more remarkable that in the full theory 
all (i dependence precisely cancels. We note that this is 
similar to what is expected to occur in phase-quenched 
QCD when m w /2 < fi < m^v/3 [2q . Beyond the critical 
chemical potential, the density increases rapidly in the 
phase-quenched theory as well. 

We have also computed the average phase factor in 
the phase-quenched theory using real Langevin dynam- 
ics. The results are shown in Fig. |U For large chemical 
potential the phase factor goes to zero on all lattices, 
making reweighting impossible. For small chemical po- 
tential, the phase factor goes to zero exponentially fast 
as the four- volume increases. This is precisely how the 
average phase factor is expected to behave [13] • However, 
we stress again that in the simulations of the full theory 
no negative impact of the sign problem was found. 

Summary and outlook - We have applied stochastic 
quantization and complex Langevin dynamics to study 
the Silver Blaze problem in the relativistic Bose gas at 
finite chemical potential, both in the full and the phase 
quenched theory. We found precise agreement with theo- 
retical expectations, and no obstacles related to the sign 
problem or in taking the thermodynamic limit. These 
results clearly stimulate the application of this approach 
to more complicated theories with a sign problem, in par- 
ticular QCD at nonzero baryon density. 

In this context two aspects have to be mentioned. 
First, due to the complexification, the Langevin dynam- 
ics no longer takes place in SU(3) but instead in SL(3, 
C). This has been discussed in detail in Ref. [Hf , where 
first numerical results for QCD with a nonzero density 
of static quarks can be found. Second, the inclusion 
of dynamical fermions, not discussed explicitly in Ref. 
[23| . is relatively straightforward. After integrating out 
the fermion fields, the fermion determinant contributes 
to the force term for the gauge fields. This is not differ- 
ent from any conventional algorithm used in lattice QCD 
simulations, except that the force is now complex, mak- 
ing necessary the extension from SU(3) to SL(3,C). For 
details concerning the inclusion of fermions in Langevin 
dynamics, see e.g. Ref. [281 ] . Work in this direction is 
currently in progress. 

Discussions with I.O. Stamatescu and S. Hands are 
greatly appreciated. This work is supported by STFC. 
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